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Abstract. We investigate in more depth the issue of backreaction in models that attempt 
at generating cosmological magnetic fields at inflation. By choosing different, physically 
motivated, parametrisations, we are able to isolate the heart of the problem, namely the 
existence, alongside the wanted magnetic field, of its electric counterpart, which turns out 
quite generally to be stronger and redder. We were also able to identify a few more interwoven 
weak spots (the typically very high scale of inflation, the width of the spectrum of modes 
processed by inflation, the blindness of the amplification mechanism to the energy scale 
processed), in a way independent on the specifications of the coupling between inflation and 
electromagnetism. Despite having stripped down the problem to the core, the obstacles 
encountered appear insurmountable, thereby posing a challenge to inflation as the incubator 
of cosmological magnetism. 
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1 Introduction 

The origin of large scale cosmic magnetic fields pertains to the collection of unsolved problems 
of cosmological particle physics. Such fields manifest themselves in a variety of environments, 
from galaxies to clusters to filaments and beyond, with very similar intensity, usually ranging 
from a few uGauss to hundreds of /uGauss; they have surprisingly been detected also in regions 
of the sky, the intergalactic medium, where little to almost no gas is present. Some useful and 
comprehensive reviews can be found in [1-15], where a detailed account of the multifaceted 
aspects of their detection and of the attempts made toward an organic explanation of their 
birth and development is given. 

While several, more or less contrived, more or less inventive, models relying on purely 
astrophysical phenomena and resting on conventional field theory can be built with various 
degrees of success, and while phase transitions constructed around standard or non-standard 
particle physics theories do give rise to magnetic seeds with very short coherence lengths 
to be fed to the magnetohydrodynamical plasma, mechanisms operating during an early 
inflationary epoch of expansion by means of a, direct or indirect, coupling between the field 
responsible for the de Sitter evolution and photons remain a theorist's favourite, due to 
relative mathematical simplicity, the obvious appeal of inflation as an anyhow necessary 
mechanism to generate primordial perturbations, and a natural, automatic solution of the 
correlation length issue. 

Indeed, in surveying the copious variety of models and constructions available, what 
catches the eye of a reader is the sheer difficulty encountered in trying to stretch high inten- 
sity magnetic fields from short to large wavelengths, either through plasma power redistribu- 
tion mechanisms (inverse cascade) or through violent astrophysical phenomena (ejection of 
magnetic field lines in the outer space) - see [16] for a possible exception. Infiation, on the 
other hand, has it all from the beginning, because any quantum fiuctuation will be pulled 
apart (causally) outside the Hubble scale during de Sitter acceleration, and will comfortably 
re-enter during a later slowing expansion stage of the evolution of the Universe. 

Generating an electromagnetic field during inflation needs however some beyond the 
standard model input in order to operate, because electromagnetism in four dimensions is 
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conformally invariant, and our Friedmann-Lemaitre- Robertson- Walker description of the Uni- 
verse is conformally flat: the combination of the two results in conformal triviality, meaning 
that the Euler-Lagrange equations of motion for the electromagnetic vector potential can be 
cast in their Minkowskian form upon a conformal transformation, and would never notice 
the presence of the non-trivial, time-dependent background. The only exception to this is 
the quantum conformal anomaly [17], which unfortunately does not appear to be viable. 

In spite of all these encouraging news, inflation as a locus for the dawn of magnetic 
fields suffers from one major drawback, that of backreaction. What typically happens is 
that, although initially negligible, the effects of the growing electromagnetic waves onto the 
background evolution turn out to be so important that de Sitter expansion could not be 
supported for a period of time long enough to address successfully the problems inflation was 
brought in to solve in the first place. What one finds is that in a large class of well justified 
models, either the produced magnetic fields will be strong, but inflation would end abruptly 
early, or inflation survives the threat but the magnetic outcome is miserably weak to be of 
any help in reconciling with observations [18-20] ^. 

In this work we want to dwell deeper on this apparently insurmountable issue, in an 
attempt to better understand and isolate the details which ultimately cause the wreckage 
of this approach; we wish to proceed in a model-independent way by parametrising the 
growth of the electromagnetic vector potential without relying on the specifications of the 
background particle theoretical description; in light of the physics thereby unveiled one may 
hope to recognise possible and potential avenues circumventing the obstacles obstructing the 
way to inflationary magnetic fields. 

The procedure we adopt in this paper is as follows. First, in Sec. 2 we will describe the 
type of theories we deal with, and swiftly skim over the arguments which show the rise of 
backreaction; from there we move on to the different parametrisations we study, correspond- 
ing to different underlying physics, in Sec. 2.1, and describe the background dynamics in 
Sec. 2.2. With all the ingredients at hand we show how much magnetic field can possibly be 
generated in a variety of situations in Sec. 3, detailing on different perspectives in tackling the 
subject in Sees. 3.1 and 3.2; the discussion of these results is given in Sec. 3.3, while an inter- 
esting promising, but ultimately doomed, counterexample is outlined in Sec. 3.4. We review 
and summarise the findings of this work in Sec. 4, alongside proposed future developments. 

2 The development of electric and magnetic fields 

Generally inflation and electromagnetism are coupled by means of one of the three interaction 
terms 

^^77 = + \f{v)FF - \rn{^?A' , (2-1) 

where the F^i, = dfj^A^, — d^A^ is the conventional electromagnetic field strength tensor, 
pfiu ^ i^fiuXnp^^ (with r^^^T^ being the totally antisymmetric Levi-Civita tensor with 7^0123 _ 

Y^l — g\) is its dual - we adopt mostly negative metric signature (+,—,—,—). The non trivial 
dynamics comes from the time-dependence of the functions I{ri), f{r]), and ni{rj) where rj 
denotes conformal time and the background inflaton field ip is hence a function of time only. 
The equations of motion for the electromagnetic field, which look like 

V^Fi: + 2{d^ \nI)FI} - {d^f )FI^ + m'A, = , (2.2) 
'^Backreactions were first mentioned in an earlier work by Campanelli [21]. 
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encode the modified dynamics of tlie system. In general terms one can then proceed to 
parametrising the evolution of the I^, /, and m?, as a function of conformal time, or the 
scale factor, and investigate the kind of behaviour the electromagnetic four-potential in 
a given gauge. The simplest way is to quantise the field (we adhere to the definitions of [22]) 
and look at the specific evolution, for instance, of the helicity /i = it Fourier modes Ah - we 
will always work in Coulomb-transverse gauge for which (A^) = {0,Ai) with diAi = 0. 

We instead opt for a different approach, since we are ultimately interested in how 
strong a generated magnetic field could possibly be, where the only constraint we have is 
that inflation develops and completes successfully. We therefore wish to be able to analyse the 
evolution of the Fourier modes Ah detachedly from the details of the form of the interaction 
term; we seek to bring to the forefront the places where the backreaction comes in, and 
eviscerate the physics at the core of the problem. In particular, we want to understand 
whether there are, in principle, ways for the Ah modes to evolve developing a magnetic field 
at large scales which is appreciable today; we aim at clarifying the relative weights of the 
electric and magnetic fields, and that of the different modes involved; we further look at the 
role of the magnetic (and electric) power spectrum in determining the final result, including 
the possibility of explicit cutoffs beyond those traditionally assumed in coping with this issue. 

Before moving on to the actual discussion, the key definitions we will employ in the rest 
of the paper are the magnetic power spectrum 

^B = 2^^^^ (2-3) 

h 

where k is the wave number of the given Fourier mode Ah, and o(r/) is the scale factor of the 
Universe. The energy density in electromagnetic field is expressed as 

h 

where the first term is related to the electric field, and the latter is the magnetic energy 
density proper. The vacuum normalisation for Ah is 

AT" = ^e-^'^" . (2.5) 



During the quasi de Sitter expansion epoch, conformal time, scale factor, number of 
e-foldings and the various Hubble parameters are related as 

a = aie^ , = rue'^ , N = ! Hjdt , 'H = - = aHi = a, (2.6) 

J a 

where t is cosmic time, the scale factors aj, a^, and tell us about the beginning and end 
of inflation, and the moment at which a given mode exits the Hubble scale = k/Hi, Hj 
being the approximately constant Hubble parameter, which we normally compute at the end 
of inflation. Also, recall that during inflation conformal time spans the range [— oo,0], and 
Ti = —l/rj. Finally, conformal time derivative is denoted with a prime, the dot being reserved 
for cosmic time. 

Already at this stage we can obtain a very primitive estimate for the suppression of the 
final magnetic strength in the most basic case. Indeed, without specifying anything about 
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the coupling in (2.1), and the details of the evolution of Ah, we find that the most favourable 
scenario only yields around 10"^*^ Gauss on Mpc scales today. This can be seen as follows. 
Imagine the spectrum (2.3) is some power of k times an amplification factor due to the 
interaction with the inflaton: 



k 



4+n 



(2.7) 



which translates in magnetic and electric energy densities 



p% ~ Amp 
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respectively; notice how the electric part is either dominating or equal to the magnetic one. 

There are two distinct possibilities, for either n > —2 or n < —2. In the first case the 
ratio between the photons energy density and the background one is p%l p^p — Amp^iJ|/M|, 
which has to stay below one. This implies a maximal amplification factor of Amp = M^/Hi, 
and, assuming radiation dominated the energy budget of the Universe at Oe, a magnetic field 
strength today 



_k_ 



4+n 



— Pcmb 



k 
He 



A+n 



(2.9) 



If we translate this in Gauss units we find approximately 

^ \2+n/2 



51^ 



10" 



Gauss , 



(2.10) 



where the most favourable case is for n — >• — 2. The very wide stretch between large scale fc, 
say 1/Mpc, and He then translates in the 10"^*^ Gauss upper limit mentioned before. The E 
and B spectra in this case are both very blue. 

The other possibility is that n < — 2, which means that the maximal amplification factor 
is Amp = (ae/ai)^+"M4/ff/; hence 
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(2.11) 



where in the last step we have reasoned that inflation has to last enough to at least include 
the scale 1/Hq. Again, in numbers this means 



6% ~ 10-"^ 



k 

He 



Gauss . 



(2.12) 



To minimise the suppression factor k/Ho the best choice is once more n — )• —2, for which (2.9) 
is recovered, with identical outcome. It is important to realise that the magnetic spectrum 
ultimately wants to be blue (n = —2) which is why large scale fields are weak; this constraint 
comes from the rampant electric field, as we will discuss in details in Sec. 3.3. 
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2.1 Parametrising the evolution of A 

Model 1. Complying with the statement of intentions enunciated above, we identify four 
parametrisations which capture different aspects of the problem. The first possibility is the 
most intuitive one, and includes the analysis of [18-20]: 



Ah 



1 



l9, 



out 



-ikrj 



(2.13) 



'2k l\ai/ V«j y V«fc, 

where we introduced the Heaviside functions iJin = '&{\krj\ — l) and T?out = i}{l — \kr]\) to select 
the modes that are inside or outside the Hubble scale, respectively. The parameters set {p, q) 
describe the growth, or decay, of the vector potential in these two cases; this parametrisation 
is therefore mode-blind, in the sense that all modes are treated equally, except for their 
time spent on the two sides of the Hubble horizon, which is obviously fc-dependent. The 
{p, q) = (0, q) case is the most common one arising from interaction terms of the first and 
third type in (2.1), which normally is inefficient for subhorizon modes; here we have extended 
it to include non-trivial evolution within the Hubble scale as well. A convenient way to 
rewrite (2.13) is 
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-ikrj 



\p-q^pN ^ 



(2.14) 



where we can further rework the last term as \kr]i\P~''e'^ = \kr]Y 

At the end of infiation the power spectrum (2.3) and the total electromagnetic energy 
density (2.4) will be given by 
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k 

where Nt determines the total number of e-foldings of infiation and Tie 
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(2.15) 



(2.16) 



where the magnetic field energy density is half of the first term in the curly brackets. Notice 
the different exponents in the two terms: they will turn out to be crucial in determining 
what backreacts and when. The limits of the k integration are the maximal and minimal k 
accessible by inflation, that is, /cmax = CLeHj = Tie and /cmin = Oj-ff/ = "Hi- 

Model 2. The second possibility is related to the axial-type coupling fFF, which is typi- 
cally efficient only for a short window of time around horizon crossing due to the structure 
of the mode equation, which includes a friction term of the form kf'A!^ [23]. We choose to 
parametrise this model by 
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(2.17) 



where the Heaviside and window functions are 



l?p,e = 'd{\kr)\ - i 
W = ^pre^post = 
9 = 1?(apre 



■pre), '&post = -^iCpost - \kr]\) 

= ■!?(cpre - \kr]\)'d{\kr]\ - Cpost) 
-ai) , @ = ^{ai - apre) 



(2.18) 
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We use the parameters Cpre > 1 and Cpost < 1 and the relative Heaviside unit steps 'dpre and 
^?post to describe modes which are not close enough to the Hubble scale, or are far enough 
out, to be not affected by the growth parametrised by r, and confined to the window W. The 
additional Q and @ are needed in order to ensure that the amplification window does not 
begin before ai. We have further defined a^re = k/cp^eHj and Opost = Cpostk/Hf. In what 
follows, for simplicity, we will choose Cpre = c = l/cpost! notice that we demand that at least 
once the full window Opre to Opost is traversed by inflation, which implies that Oe/c > coj, or, 
equivalently, that < e^*. Hence, the parameters set we work with is (r, c). 



Again, it proves convenient to rewrite (2.17) in terms of N as 
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(2.19) 



from which we can compute the power spectrum as 
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(2.20) 



where the first term refers to high energy modes for which the amplification window was 
not traversed completely due to the end of inflation, the second term encompasses all modes 
which have been fully amplified, and the last, third, term, describes low-energy modes which 
were being amplified right at the beginning of inflation. The final energy density is given by 
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(2.21) 



where the magnetic field energy density is made of half the first, second, and last term, 
the rest pertaining to the electric field. We can understand the first piece as belonging to 
the highest energy scales, which were not fully amplified, the second piece represents all 
intermediate scales with full amplification, the third is similar to the first but arises from the 
conformal time derivative in the electric field, and the last one accounts for the lowest energy 
modes which were also not fully amplified. Notice that the result is barely dependent on the 
total number of e-foldings as it should be, as the width of the window is independent on the 
dynamics of inflation. 

Model 3. The parametrisations described so far are not knowledgeable of the mode being 
processed, except for its being in, out, or around the time of horizon crossing. This is 
normally the case as the coupling functions are independent of k itself. However, one can 
envisage a situation where this restriction is dropped, and the efficiency at which a mode is 
amplified does depend directly on the mode itself. Since the typical amplification mechanism 
is operative outside the Hubble scale, we choose to allow for a /c-dependence in the exponent 
q, which becomes now q = q{k/k)'^, and where the reference scale k is fixed at 1/Mpc, for 
this is the characteristic scale we will explore numerically. 
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The parameters set thus becomes {q,d), with d brought in to define the "tilt" of the 
amplification mechanism; positive values for d mean that the higher end of the spectrum will 
have a higher boost factor compared to the red end, while a negative d would favour large 
scales over blue ones. Explicitly: 




(2.22) 



or, in terms of e-foldings, 

A = -i= I'&in + \kv\~^'&ont} e" 
v2k ^ 

Once again, the final power spectrum at the end of inflation is 

and the total electric and magnetic energy density is the integral 

= ^ y" dK {2^3 + nf} , (2.25) 

where we have defined k = k/Tie for convenience. It is in general not possible to solve (2.25) 
in a closed form, due to the /c-dependent exponent q, so we will limit ourselves to directly 
solve the integral by numerical means. 

Model 4. This parametrisation deals with the issue of explicit cutoffs in the magnetic (and 
electric) field spectrum. So far one further, unspoken, but basic assumption has been that 
all modes from k^in up to /cmax partake in the amplification process. As is clear already by 
peeking at the simplest results Eqs. (2.16) and (2.21), the highest scales are very likely to 
be dominant in the energy budget of the electromagnetic field, and the most dangerous in 
terms of their possible backreactions. The fact that the total energy density has to be shared 
by such a wide range of modes makes it much harder for single modes, especially large scale 
modes, to be entitled to any relevant portion of it. One may thence hope that by cutting off 
part of the spectrum, because of some physical cutoff in the theory which produces the 
interaction term, the situation could be ameliorated. 

The range of modes which appear to be the trickiest to generate has been already 
identified in the low energy one, at and around the Mpc scale. Thus, of the two possible 
cutoffs, the infrared and the ultraviolet one, the most interesting would be the latter, for there 
are only a handful of modes below 1 /Mpc which could be removed. Physically speaking, it 
also seems reasonable to work with an ultraviolet cutoff for its connection to the effective field 
theory producing the higher order coupling terms in the Lagrangian. While we will retain 
the p = assumption in this class of models, we will discuss both cutoffs in what follows, in 
connection to their different physical interpretation. 

The high energy cutoff in the theory is implemented by simply restricting the range of 
integration in (2.4), and attaching a Heaviside step to the power spectrum (2.3) as 
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(2.27) 



where we have used < c < 1 to set the upper hmit for k through i^uv = '^{cl-Le — k). 
The parameters set is ((7,c), however, these expressions are easily generahsable to include 
within- horizon amplification (through p) , or accommodate an infrared cutoff (through = 
'&{ck — T-ii)): in the latter case the final energy density (2.27) is modified as 
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(2.28) 



Model e. In this last example we want to see what happens in the extreme case of an 
exponential growth for Ah- The relevance of this possibility lies in the fact that it is perhaps 
the most immediate and direct way to put electric and magnetic field on the same level, if 
not even to suppress the former compared to the latter. The details of the struggle for power 
of the two will be analysed below, so for now let us simply introduce the parameters {p, q) 
or (P, Q) = {puj, quj) as 
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(2.29) 



where we can think of the frequency co as exactly k or some independent parameter (thence 
the re-parametrisation pu ^ P and quj ^ Q). 

For the last time, let us also explicitly write down the final power spectrum (2.3) 
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the total electromagnetic energy density (2.4) reads 
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(2.30) 



(2.31) 



(2.32) 



The electric field contribution here, besides the usual first half of the first term in the curly 
brackets of (2.32) which matches that of the magnetic field, is proportional to the square of 
quj/He, and will depend on the specific values for lo and q. This expression can be simplified 
by introducing Q = Q/He, P = P/He, and ignoring the \rjk\ <C \r]i\ and |r/e| <C \rik\ in the 
exponents, to give 
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(2.33) 



The last expression is technically valid for values of Q = quj > Hi, but due to the extra 
factor of exp(— 4A'^t) in the second term, whenever Q ^Hi the last term becomes completely 
insignificant anyhow. 
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2.2 Background dynamics 

The background evolution is determined by the details of the inflationary expansion epoeh, 
which we characterise in this section. Although we will stick to a very simple and circumstan- 
tiated example (mostly chaotic inflation), the results we obtain do apply generally to other 
types of inflation; in the end the only parameter which makes the difference is the scale of in- 
flation, which sets the "distance" between the end of inflation to today, thereby determining 
the dilution factor for the magnetic energy density and strength at given wavelength. 

The background dynamics is governed by a scalar field the inflaton, with potential 
V{ip). In order to be definite, and in order to be able to fix all parameters for the subsequent 
numerical calculations, we choose to work with either a monomial power-law potential or the 
type 

V{ip) = (^j^y , for a 7^ 4 or V{ip) = A(^^ for a = 4, (2.34) 

where the mass scale M (or the coupling constant A) determines the scale of inflation; we 
will also discuss briefly one example of small-field inflation, exemplified by the potential 

V{(p) = M"^ {1 + cos gip) , (2.35) 

where the dimensionful coupling g is the inverse "decay constant" for a would-be (pseudo) 
Nambu-Goldstone boson {(p in this case). 

Given the potential one can readily obtain the dynamics, assuming there are no back- 
reactions in effect: 

^ + SHif + V^ = 0, (2.36) 

where = dV/d(p. The inflaton will start somewhere up high in the potential, and slowly 
roll down until its kinetic energy density, which is assumed to be unimportant at this stage, 
catches up with the decaying vacuum potential energy, and finally causes the inflaton to 
decay while rapidly oscillating around its minimum. A handy way to follow this evolution is 
through the number of e-foldings for which inflation lasts, which we can approximate as 

iV=^/^d,. (2.37) 
The end of inflation is typically estimated to happen when the slow-roll parameters 

_ H Mi fV^Y _ isR ,.2fy<pA ^o',^^ 



become of order one, hence ^Pf ^ aM4/y/2 (similar results hold for the a = 4 case). This 
means that the evolution of itself is 

(^f«^^+2«iV. (2.39) 

The scale of inflation is determined by M, which ultimately will yield Hj, wherefrom 
we can make the connection to our Universe by assuming the type of evolution (matter 

radiation, or else) in the wake of the aftermath of the de Sitter expansion. This parameter 
is fixed, for a given Nf by normalising the generated curvature perturbations to the directly 
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related temperature fluctuations as measured by CMB experiments such as WMAP. The 
power spectrum of such temperature fluctuations is encoded in the quantity 

6j, = \ , « (2 X W-'f , (2.40) 

where the field (p is evaluated at the end of inflation. Using this result, one can fix the scale 
of inflation H] = V/3M^ as 

H] = ^2^!^^4m| . (2.41) 

^ (1 + 4iVj/a)°/2+i " ^ ^ ^ 

One more piece of information which we will be using, effectively as a way to parametrise 
variations in the total life span of inflation, is the spectral index of the aforementioned 
perturbations, which for monomial inflation reads 

„. = l_Mei±|. (2.42) 

This slightly non-standard take on inflation is adopted for the ease of numerical im- 
plementation which thereby comes about; the results we will obtain and document do not 
crucially depend on these details however, and can be straightforwardly adapted and gener- 
alised. In our analysis we will always work with the total Nt also because we do not want 
to descend too deeply in the swampy lands of actual quantisation and renormalisation of the 
electromagnetic potential (and in turn of the inflaton field [24, 25]), which will take us far 
beyond the scope of this paper. We will hence be content with always "booking" a range 
of modes going from /cmin to fcmax which depends on Nt in its width, and on Hf (which can 
or can not be constrained given Nt) for its absolute value; clearly, one needs to ensure, a 
posteriori, that the modes we are looking at are within the range touched by inflation. 

3 How strong can the magnetic fields possibly be? 

The key quantity which roughly determines the onset of backreactions is the ratio 

x=^; (3.1) 

when this ratio becomes of order one (being negligible at the beginning of inflation) the 
background becomes unstable, since the Friedmann equation — + /O7 does not admit 
quasi de Sitter solutions anymore. Obviously we do not really know what happens when 
X — 7- 1 until we solve the coupled Euler-Lagrange equations for ip and Ah, but it seems safe to 
assume that, if inflation surpasses the impediment, the final magnetic field will be negligible, 
see [18-20, 23] 2. 

The time at which the Ah backreacts could be found instead by directly looking at the 
equations of motion (2.36) where terms explicitly containing Iip{i]), fip{f]), or m(r/) appear 
multiplying some combination of Ah- The balance is then between V^p and the latter, and can 
in principle be very different from x. However, to investigate this we need to know explicitly 
the form of the interaction term, so we do not pursue this possibility any further, bearing in 
mind that our results are a conservative estimate, in the sense that the upper limits we will 
derive could be further lowered. 



^Notice that in the specific model of [21] the problem does not arise because only large scales are touched, 
and the amplification takes place while these are superhorizon: as we find here, this might be the only viable 
combination. 
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3.1 Finding the maximal allowed strength 

Let us pick a particular background and look at the limits on the model parameters. To 
begin with, we will choose to fix the spectral index of (2.42) to Us = 0.97 and normalise the 
density fluctuations 6h as in (2.40). This for a given potential, say chaotic inflation with 
a quadratic potential a = 2, provides us with the total duration of inflation, in this case 
around Nt « 66. The only unknowns are then the model parameters sets: {p, q), (r, c), {q, d), 
and {q, c), for the four models we have introduced above. 

We look at regions in the parameters space where the ratio (3.1) stays below one, 
which are plotted in Figs. 1, 2, 3, and 5; on top of the allowed and disallowed regions 
(according to this criterion) we draw the contours representing (the logarithm of) the strength 
of the magnetic held 6^ today. Notice that, again for the sake of simplicity, we assume that 
the Universe reheats perfectly and instantaneously, delivering a beautiful radiation-flUed 
primordial soup; this is expressible as 

p'^^c^^T^, ae = ^ao, (3.2) 

e 

where oq = 1, and Tq = 2.75 K is today's cosmic microwave background's temperature - we 
are disregarding the contribution of neutrinos, as well as the latter dark energy. 
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Figure 1. Allowed (in light blue) and disallowed (in light red) regions in the {p, q) parameters space. 
Superimposed are contours of constant log 5^ today, at 1/Mpc (in Gauss units). 

Model 1. The first parametrisation, Eq. (2.14) plainly describes a simple growth inside 
and outside the horizon at different rates p and g, respectively. The allowed regions and 
the generated magnetic field strength are shown in Fig. 1. Clearly smaller values of these 
parameters will be required for a long enough inflation. Notice that the parameter q can never 
be bigger than approximately 1.2, meaning that the magnetic held spectrum will always be 
blue. 
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Figure 2. Allowed (in light blue) and disallowed (in light red) regions in the (r, logc) parameters 
space. Superimposed are contours of constant log (5^ today, at 1/Mpc (in Gauss units). 

Model 2. The second parameters set, (r, c) from (2.19), accounts for a fixed-time growth 
around horizon crossing (as it typical in models where photons and inflation are axially 
coupled). As expected, see Fig. 2, large r are not allowed, especially as we widen the window 
with c, since in that case the modes would be amplified for longer, and backreact more easily. 
We have also noticed already that this result depends only very little on the total number 
of e-foldings. The spectrum is typically blue, although it is possible to have a red "tail" for 
high energy modes: the large scale modes are in any case very blue. 

Model 3. The running of the amplification factor with k is expressed in (2.23) through the 
set {q,d). Notice (Fig. 3) that the largest magnetic field strength on large scales is obtained 
for d ~ 0, meaning that uniform amplification is preferred. We can understand this by 
recalling that the range of modes in question is numerically huge, as k E 1], and even 

relatively small values for d will enhance dramatically one or the other end of the spectrum, 
thereby tightening the constraints on the growth parameter q. Looking at it in a logarithmic 
scale in Fig. 4 we see that very moderate values for d in fact do ameliorate the situation, but 
clearly not enough to allow for a significant magnetic field today. 

Model 4. The two following examples refer to Eqs. (2.27), and (2.28), Fig 5, left and right 
panels, for the ultraviolet and infrared cutoffs, respectively. In the first case the change 
is minimal, because once again, due to the different powers for the electric and magnetic 
contributions, as soon as g > 1 the electric spectrum turns red, and the high energy cutoff 
becomes uninfluential. In the opposite case instead, where we cut off large scale modes, there 
can be seen a fairly more significant difference; nonetheless, this possibility is still not viable, 
because there are only a handful of modes which can be cut off before reaching 1 Mpc, and 
at this scale the electric field is still too strong, and redder, than the magnetic field. 
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Figure 3. Allowed (in light blue) and disallowed (in light red) regions in the (d, q) parameters space. 
Superimposed are contours of constant log (5^ today, at 1/Mpc (in Gauss units). 

Model e. Speaking of the relevance of the electric field, as will be sketched in the next 
section, one simple way to prevent it from growing much faster than the magnetic field is to 
engineer some exponential-type of amplification for Ah', this is what was proposed in (2.29). 
The spectrum (2.31) is blue unless the exponential tilts it; in particular, the lowest energy 
mode /cmin = has the same power of the highest fcmax = T~{-e (or more) if 2q(jj/kyam > ^Nt, 
even though it decays exponentially for higher modes, until the overcomes it. A closer 
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Figure 4. Zoom on Fig. 3 for small negative (left panel) and positive (right panel) values of d, 
expressed in terms of — log(— d) and log(d), respectively. 
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look at Eq. (2.33) reveals that as soon as the spectrum "turns red", the electric contribution 
becomes dominant once again, because of the qui/T-Li prefactor, which is large in this case: 
the electric field rears its ugly head again. Despite the ratio x being still well below one at 
the turning point for quo, it climbs up exponentially rapidly, and does not allow for more than 
about 10~^^ Gauss today on scales of a Mpc. 

3.2 Variations on a theme 

There are several variations on this work which can be investigated; among the several 
possibilities available we focus on a few examples. 

Varying n^. First, we can vary the scalar index Ug, which effectively is a way to play 
with the total number of e-foldings. Of course, as long as we stick to the simple quadratic 
potential, and given that we need at least around 60 e-foldings to be able to let today's 
Hubble scale T-Lq be inside the inflating machine once, there is only that much freedom in 
choosing Ug. In Figs. 6, and 7 we compare the allowed regions where a; < 1 where we rise 
to 0.98 and 0.99, respectively. We do so for the first two models with sets {p,q) and (r, c), 
respectively; in the first case the "good" parameters space further shrinks as expected - more 
inflation means more growth; in the last case the results are barely dependent on Nt, if not 
for a slight widening of the blue region. 

Varying V{(p). Next, we can tamper with the inflaton potential with the hope of finding 
different relations between the background parameters n^, Hj, and Nt, which could help in 
alleviating the severe dilution the generated magnetic field incurs in. Since we have not found 
any significant departure from the already reported outcomes by simply switching to quartic 
or other simple monomial potential, we turn our attention to the "natural" inflation case, 
Eq. (2.35). 
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Figure 5. Allowed (in light blue) and disallowed (in light red) regions in the (log c, q) parameters 
space for the ultraviolet (left panel) and infrared (right panel) cutoffs. Superimposed are contours of 
constant log 5^ today, at 1/Mpc (in Gauss units). The dashed line in the right panel marks the value 
of c for which 1 Mpc is cut off; no power on this scale is generated for values of c below this line. 
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Figure 6. Same as Fig. 1, but for Ug = 0.98 (left panel) and rig = 0.99 (right panel); this changes 
the overall duration of inflation to Nt — 99.5 and Nt — 199.5, respectively. 



In natural inflation the scalar index turns out to be directly linked to the inflaton's 
inverse decay constant g by 

= 1 - (gM^f , (3.3) 

which means that if observations pin down the value for Ug, we automatically fix g, but we 
still have complete freedom as for Nt, unlike the case of a monomial potential, see (2.42). The 
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Figure 7. Same as Fig. 2, but for Ug — 0.98 (left panel) and Ug — 0.99 (right panel); this changes 
the overall duration of inflation to Nt = 99.5 and N — 199.5, respectively. 
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most important consequence of this is reflected in the final Hj which now explicitly depends 
on Us and Nt separately: 

^ ^_^s2^^{ns-l)N,Mi , (3.4) 

up to a factor of order unity, weakly dependent on Ug itself. However, with typical values 
for the scalar index and the total e-foldings, also (3.4) falls in the usual range of around and 
above 10^^ GeV, and no significant changes occur. 

Flipping the argument. One can embrace the somewhat backwards approach of deciding 
first what the desired magnetic field strength is, and given that compute the allowed duration 
of inflation for given set of parameters. For instance, we can agree to succumb to the 
observational imperative demanding nGauss strength on scales of a Mpc, and work out 
how long a successful inflation could possibly last. This point of view is in fact a simple 
"re- normalisation" of the results of the previous sections, for instead of asking how large will 
6b he with the constraint x = 1, we ask how long will x < 1 be satisfied for given the final 
6b, and we will not pursue it further. 

3.3 Some physics 

The results of the previous section show how undoubtedly challenging it is to generate mag- 
netic fields of tangible intensity today, the more so for large scales, which is where inflation 
is expected to be the most helpful. We had undertaken this investigation with the hope to 
shed some light unto the details for which this is so, and to be able to point out which of the 
gears operating in the mechanism are responsible for its ultimate failure. 

The role of E. First, in all cases analysed above we have observed how the different 
powers appearing in the expressions for the electric and magnetic energy densities result in 
the prevalence of the former over the latter, and its ultimate fundamental role in stopping 
inflation while 6b is still small. We can schematically formalise this finding as follows. The 
electric and magnetic fields, in Fourier space, are related to the vector potential Ah as 



\dr,Ah\\ Bl:^k^\Ah\^. (3.5) 



Normally we parametrise the vector potential as Ah = viv)^ ^^^^ which means that 



1 + 



ky 



Bl . (3.6) 



Now, typically y' ~ y/?7, which translates in y' /ky ~ ^/kr], and once plugged into the 
expression for p^, Eq. (2.4), says 



Pe- J dlnkEl'^ ~ J dlnfc 



1 + 



kv 



k 



Bl'^ . (3.7) 



We can formally write f din kB'^'^ = B^'^. Normally we also parametrise the fc-dependence 
of Bk in terms of a power law, for instance A;", and, up to some order one coefficients, we 
obtain 

p% B^Z. - + C{B^^l. - e^'^'B'jJ 

c^il- e-2^^0i?^ax + C{1 - e2(i-^)^0i3^ax , (3.8) 
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where the 'max' and 'min' subscripts as usual refer to the maximal and minimal mode to 
be taken into consideration and C is an order 1 constant of integration; notice that at the 
other end of the spectrum B'^^^ = e~^^^ B'^^-^. The first term is common to both the electric 
and the magnetic field, whereas the second piece is peculiar to pe only. If we pick a blue 
spectrum for the magnetic field, that is, 6 > 0, we see that either the electric field is also blue 
(for h > 1) and their energy densities are comparable, or, at smaller h it is red, and steals 
most of the energy density of because of the last exponential. If the magnetic field has a 
red spectrum then, which is what one wishes to fit observations, the p^^ is always dominating 
by a factor of e^^* . This is the source of the problem: one can not simply get away with some 
sizeable, red tilted, magnetic field energy density without a much stronger electric field. 

The possible way to avoid this is to consider coupling constants which result in exponen- 
tial amplification of the Ah, leading to y(r/) ~ e'^^ or some more complicated exponentials, 
in which case the precise form and coefficients appearing in the exponential will determine 
the winner in the electric versus magnetic competition. For instance, a Gaussian type of 
evolution, y ~ exp(fc7/)^ would reverse the situation in favour of the magnetic field. 

Low scale inflation. One possible place of intervention to cure these models would be 
the overall scale of infiation, parametrised by either Hi or Og. Indeed, a much lower scale of 
inflation would imply much less dilution until today for the magnetic power, which in turn 
translates into less magnetic field wanted at the end of inflation. 

However, as we lower the scale of inflation, while the required magnetic strength (squared) 
(5^ decreases as l/Og, at the same time the flnal energy density of inflation also decays. 
As long as the Universe was always dominated by radiation after inflation, and the reheating 
process is instantaneous and has no effects on the background dynamics, then also p^ follows 
the l/Cg law, thereby nullifying the positive effects of the lower Hi . Notice that the late-time 
matter dominated era does not change significantly this result, for which instead a signif- 
icantly long lasting period of, say, matter dominated reheating (or any other background 
behaviour slower than 1/a^) would be needed. 

One further slight improvement brought about by this possibility is that the required 
minimal number of e-foldings lowers, since there is less of a stretch between the highest 
(now lower) mode and today's Hubble scale. As it is clear from the figures of the previous 
sections however, such small, in terms of the actual values involved, adjustments do not 
appear feasible in restoring the effectiveness of these models - we are attempting at bridging 
many orders of magnitude gap in (5^ by means of order one or ten corrections. This avenue 
hence does not seem to be harbinger of any significant progress. 

3.4 A spiky counter-example 

In light of the insight thus far developed, one is led to believe that if we were able to cut 
the spectrum of amplified modes down to a very narrow band centred around the modes of 
most interest, then we could be able to work our way around backreactions; moreover, one 
would like to keep the electric field under control while still allowing for a large magnetic 
field. These two characteristics are naturally obtained from models where the growth of the 
Ah is resonant, that is, Ah — e^^^"^ , where the time span Ar/ depends on the details of the 
coupling and uj is the frequency of the resonance. Such models can be constructed from axial 
couplings quite easily [26], and have the further attractive feature that the coupling constant 
never needs to be large for the resonances to be efficient ^. 

^Interesting bounds on this type of coupling come from non-gaussianity, see, e.g., [27]. 
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Despite this being the most favourable situation one could possibly imagine, the fact 
that the resonance is confined to be operative within the Hubble scale is an inherent issue 
which ultimately will determine the demise of this kind of models. To prove this, we will 
assume that 

1. the amplification mechanism is operative only inside the horizon, and ceases, with 
consequent freezing of the vector potential, once a mode becomes superhorizon; 

2. the only effect of an eventual reheating stage following the de Sitter epoch is only to 
alter the overall decay of the total infiaton energy density, and has no impact on the 

electromagnetic field; 

3. backreactions are of the background type (they impact the Friedmann equations) - this 
unlaces the following estimates from a particular coupling term, and is a conservative 
hypothesis; 

4. the spectrum is a spike at a given wavelength; notice that this represents the best case 
scenario, since any different spectrum would necessarily have a higher energy density 
associated with it (at the same power for the pivotal wavelength) and would make the 
backreactions problem worse. 

In the case of a delta function type of spectrum, the magnetic energy density pB and 
the power spectrum (5^ coincide up to a factor of 4, and we can use either one of them in our 
estimates. Now, assume that we are interested in a given final value for (5^*^ stored primarily 
at a scale which, for illustrative purposes, we take tohe k = k = 1/Mpc. If nothing happens 
to the magnetic field from the end of inflation until today, we can estimate its value as 

The value of is what we would like to obtain from the interaction with the rolling infiaton, 
and can be chosen to be the observed value at that particular lengthscale, in this case 1 nGauss 
(squared), or a smaller value which would trigger, e.g., a dynamo amplification, say 10^^*^ 
Gauss. If this field is frozen from the moment it is stretched beyond the Hubble horizon (as 
per assumption), this field will continue its evolution backwards into the inflationary epoch 
as 

Pb{N) = p^e^"" . (3.10) 

The scale of inflation, determined by Oe, can be estimated using entropy conservation; 
we can allow for a period of matter domination behaviour after the end of inflation (reheating) 
when the scalar field is oscillating, and assuming all energy density of the Universe originates 
from the decay of the field itself, one can estimate 

e ~ ( ^^iV^ ( — ' 

~ \ar) Xao) 

where a,, refers to the scale at which radiation domination kicks in. If we then look at the 
end of the de Sitter stage, we find a ratio 



(3.9) 



(3.11) 



This ratio, going back into the actual accelerated expansion epoch, will, within our approxi- 
mations, grow as the fourth power of rj = rjeC^ . The final ingredient necessary is the number 
of e- folds at horizon exit, for a given mode k, and is given by 



There are two alternative ways to demonstrate that in this situation, the electromagnetic 
field (which is essentially the magnetic field because |5r,^/iP — /i ^ 1) is doomed 

to backreact before it can reach the desired values. One can decide on a particular value p^, 
follow it up until horizon crossing given by (3.13) and compare it to p^p at the same time: 
if this ratio is larger than one than it means that the electromagnetic field has to be larger 
than what is allowed to be at horizon exit, up until when its evolution is certain. Similarly, 
one can compute the number of e-folds at which the ratio x = x{N) = 1 and if this N turns 
out to be smaller than A^^ then the system backreacts. 

The second possibility is to choose the most optimistic option psiNk) = PipiNk), and 
follow the magnetic energy density down to today. This leads to 



where we have taken Pip{N) = p^ for simplicity. This expression (3.14) relates the largest 
possible magnetic field at a given scale k today with the duration of inflation, and, for large 
modes especially, clearly shows how any reasonable duration for inflation automatically leaves 
us with only crumbs of magnetic energy today: 66 e-foldings of chaotic (a = 2) inflation can 
give birth to a disheartening 10—51 Gauss at 1 Mpc today. 

4 Conclusions 

To conclude, we have analysed the impact of the backreactions during inflation, when the 
inflaton field itself is coupled to the electromagnetic vector potential. Thanks to a choice of 
physically motivated parametrisations, we were able to compare the allowed regions in the 
parameters spaces which lead to a controlled background, and at the same time compute the 
strength of the magnetic field today at any given interesting wavelength, without having to 
rely on a specific type of coupling, and to a certain extent, independently on the details of 
inflation. 

We primarily focussed on large scales, as they are the ones which appear to be harder to 
be generated via late-time mechanisms such as phase transitions or astrophysical processes. 
The results are somewhat discouraging, as it was essentially impossible to find a parametri- 
sation and a choice of values for the parameters which allowed any significant production 
of magnetic fields. We have elaborated on four possibilities: simple growth with different 
strength inside and outside the Hubble scale; amplification within a narrow window around 
horizon crossing; power-law mode-dependent growth; spectral cuts, both infrared and ultra- 
violet. In all these cases the strength of the magnetic fields today could never exceed 10"^*^ 
Gauss on Mpc scales. 

We have found that the chief reason behind this result is to be found in the presence of 
the electric field alongside the magnetic counterpart; the different power of momentum with 
which the former appears in the energy density makes its spectrum always redder than that 




(3.13) 




H 



(3.14) 
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of the magnetic field, for the same values of the parameters. Therefore, miless a very specific 
mechanism is able to differentiate between the two components by boosting the magnetic 
field while leaving the electric one alone, large scale power will always need to be suppressed. 
Also, the fact that normally all modes are amplified democratically, with no distinction 
among different energy scales, in conjunction with the width of the spectrum accessible by 
inflation, makes for a challenging selection of large scales to be appointed a significantly large 
power. While some of these considerations were available in the literature, scattered among 
a handful of papers, we have here tackled this issue in an organic way. 

This constitutes only a first step towards a more complete understanding of "what 
goes wrong" when we attempt to generate magnetic fields during infiation: being able to 
point out where the backreactions arise from in a model-independent way is a necessary 
starting point in building a successful model where this can be accomplished. Of course 
we have not analysed what actually happens when the electromagnetic field energy density 
approaches that of the background, and the ultimate fate of infiation depends on the details 
of the coupling through the coupled equations of motion. Furthermore, we have not dug 
in the reheating stage following inflation, where the magnetic power spectrum can change 
shape and intensity quite dramatically, as we leave this much more involved issue for future 
investigation. 

However, although we could conceive an anaesthetic combination of the elements ex- 
amined in this work which might succeed in the very end, our order of magnitude analysis 
showed how, on fairly general and solid grounds, if inflation were to survive despite the men- 
acing electromagnetic field, the resulting magnetic half would be completely negligible, the 
more so on large scales, thereby precluding the identification of such fields with the observed 
large scale magnetism. 
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